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Abstract 

Telomeres are specialized nucleoprotein structures that protect chromosomal ends from degradation. These structures 
progressively shorten during cellular division and can signal replicative senescence below a critical length. Telomere length 
is predominantly maintained by the enzyme telomerase. Significant decreases in telomere length and telomerase activity 
are associated with a host of chronic diseases; conversely their maintenance underpins the optimal function of the adaptive 
immune system. Habitual physical activity is associated with longer leukocyte telomere length; however, the precise 
mechanisms are unclear. Potential hypotheses include regulation of telomeric gene transcription and/or microRNAs 
(miRNAs). We investigated the acute exercise-induced response of telomeric genes and miRNAs in twenty-two healthy 
males (mean age = 24.1 ±1.55 years). Participants undertook 30 minutes of treadmill running at 80% of peak oxygen uptake. 
Blood samples were taken before exercise, immediately post-exercise and 60 minutes post-exercise. Total RNA from white 
blood cells was submitted to miRNA arrays and telomere extension mRNA array. Results were individually validated in white 
blood cells and sorted T cell lymphocyte subsets using quantitative real-time PCR (qPCR). Telomerase reverse transcriptase 
{TERT) mRNA (P = 0.001) and sirtuin-6 (SIRT6) (P<0.05) mRNA expression were upregulated in white blood cells after exercise. 
Fifty-six miRNAs were also differentially regulated post-exercise (FDR <0.05). In silico analysis identified four miRNAs (miR- 
186, miR-181, miR-15a and miR-96) that potentially targeted telomeric gene mRNA. The four miRNAs exhibited significant 
upregulation 60 minutes post-exercise (P<0.001). Telomeric repeat binding factor 2, interacting protein (TERF2IP) was 
identified as a potential binding target for miR-186 and miR-96 and demonstrated concomitant downregulation (P<0.01) at 
the corresponding time point. Intense cardiorespiratory exercise was sufficient to differentially regulate key telomeric genes 
and miRNAs in white blood cells. These results may provide a mechanistic insight into telomere homeostasis and improved 
immune function and physical health. 
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Introduction 

There is mounting evidence of an association between habitual 
physical activity and longer leukocyte telomere length (LTL) [1—6]. 
Telomeres are specialized nucleoprotein structures that protect the 
ends of linear chromosomes and progressively shorten with each 
round of cellular division. At a critically shortened threshold, 
genomic instability, replicative senescence and apoptosis ensue. 
Accelerated telomere shortening is associated with a host of age- 
related chronic diseases and risk factors [7- 1 5] . Adaptive immune 
cells regulate telomere length via the enzyme telomerase. 
Telomerase minimally consists of two core components; the 
catalytic subunit, telomerase reverse transcriptase (TERT) and an 
antisense RNA template (TERC). The shelterin complex, a 
dynamic conglomerate of six telomeric accessory proteins also 
plays a major role in telomere homeostasis [16]. White blood cells 
(WBCs) were chosen for telomeric gene analysis due to the high 



correlation between telomere length in these cells and those of 
other tissue types [17-20]. 

It is unknown whether the association between physical activity 
and telomere length in WBCs is due to the amelioration of 
oxidative stress and inflammation, the exercise-induced regulation 
of telomeric genes or a complex interplay between all three. 
Another possible mechanism is the differential regulation of 
microRNAs (miRNAs), which are known to respond acutely to 
physical exercise [21-24]. miRNAs are short, non-coding RNA 
molecules which post-transcriptionally regulate gene expression by 
binding to the 3 ' or 5 ' untranslated regions (UTR) of messenger 
RNA (mRNA). Despite burgeoning roles for miRNAs as potential 
mediators of exercise-induced adaptive processes, little is known 
about miRNA involvement in telomeric homeostasis. Epigenetic 
modifications play a key role in telomere length homeostasis [25- 
27], hTERT and telomerase regulation [28-30]; however, little is 
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Table 1. Physiological characteristics of the 22 male Participants. 



Characteristic 


Mean 


SD 


Age (years) 


24.0 


±7.3 


Height (cm) 


180.7 


±4.3 


Body Mass (kg) 


78.5 


±9.0 


BMI (kg/m 2 ) 


24.0 


±2.5 


Waist (cm) 


81.5 


±5.6 


Hip (cm) 


98.6 


±5.1 


Waist:hip ratio 


0.8 


±0.03 


Systolic BP (mmHg) 


130.8 


±11.7 


Diastolic BP (mmHg) 


72.0 


±8.3 


Resting HP, (b min~ 1 ) 


64.1 


±11.6 



SD (standard deviation); BMI (body mass index); BP (blood pressure); mmHg (millimetres of mercury). 
doi:1 0.1 371 /journal.pone.0092088.t001 



known about the role of miRNA-mediated regulation of telomeric 
genes. 

The aim of this study was to investigate potential mechanisms 
underpinning the positive association between physical activity 
and WBC telomere length. The specific aims were to investigate 
the acute effects of 30 minutes of intense cardiorespiratory exercise 
on the expression of genes involved in telomere regulation in 
WBCs, and to identify the exercise-induced expression patterns of 
miRNAs with potential telomeric involvement. 

Methods 

Ethics Statement 

All eligible participants read a plain language information 
statement outlining all aspects of the project in lay terminology. 
Informed consent documents explaining the purpose, potential 
risks and benefits of the project were then signed in the presence of 
a witness. The study, recruitment and consent procedures were 
approved by the Human Research Ethics Committee from 
Federation University Australia. 

Participants 

Twenty two healthy, non-smoking males (mean age = 24. 1 ± 1.5 
years) were recruited to participate in this study (Table 1). Specific 
health and medical history was obtained via a physical activity 



readiness questionnaire (PAR-Q) and general lifestyle information 
was obtained via a health and lifestyle survey. 

Physiological Measurements 

Participants were seated comfortably for 10-15 minutes prior to 
resting blood pressure and heart rate measurements. Standard, 
calibrated electronic scales were used to determine body mass and 
a standard free-standing stadiometer was used to determine height. 
Waist and hip measurements were taken at standardized sites 
using a 2 m metal anthropometry tape measure. 

Fitness Measurement 

Participants undertook a treadmill-based peak oxygen uptake 
(VC>2p,:ak) test using a Metalyser metabolic system (Cortex 
Biophysic, Leipzig, Germany). The test protocol started with a 
five minute warm up period at 10 kmh -1 after which the speed 
increased by 1 kmh 1 each minute. The incline remained at a 
constant 0% throughout the test. Breath by breath gas exchange 
and heart rate were continually monitored. The test was 
terminated when the participant indicated volitional failure and/ 
or when oxygen dynamics showed obvious levelling-off despite 
increases in work rate. 

Exercise Protocol 

A minimum of five days and no more than seven days after the 
fitness test, each participant undertook the exercise intervention. 



Start 30m in 
exercise 



Finish 30min 
exercise 



I 



30 min rest 



30min running at 

80% V0 2pMk 



60mm recovery 



Pre-exercise 
blood sample 



Post-exercise 
blood sample 



60 min post-exercise 
blood sample 



Figure 1. A schematic overview of participant blood sampling and exercise intervention. A baseline blood sample was taken 30 minutes 
before the onset of exercise. Participants then completed a 30 minute bout of treadmill running at 80% of previously determined V0 2 peak. Additional 
blood samples were taken immediately post-exercise and at 60 minutes post-exercise. 
doi:1 0.1 371 /journal.pone.0092088.g001 
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Table 2. Treadmill ramp test and exercise intervention data. 





Treadmill ramp test data 


Mean 


SD 


0 2 peak (mL kg min -1 ) 


49.3 


±4.7 


Maximum V'E (Lmin -1 ) 


125.9 


±12.4 


Maximal heart rate (b'mirT 1 ) 


178 


±8.9 


Maximum RER 


1.20 


±0.1 


Test duration (minis) 


12.5 


±1.6 


30 min exercise intervention data 


Average % of 0 2 peak during 
30 min run 


80.8 


±7.1 


Average V'E (Lmin -1 ) 


94.6 


±14.3 


Average heart rate (bmin -1 ) 


163.7 


±12.9 


Average RER 


1.0 


±0.1 



0 2 peak (highest oxygen consumption achieved in test); V'E (minute ventilation); RER (respiratory exchange ratio). 
doi:1 0.1 371 /journal.pone.0092088.t002 



Participants were asked to refrain from vigorous physical activity 
during the preceding 48 hours. A 10 ml resting blood sample was 
taken from the median cubital vein before the exercise bout. 
Participants then rested for 30 minutes before undertaking 30 
minute continuous bout of treadmill running at 80% of previously 
determined 0 2pea k- Breath by breath gas analysis was conducted 
throughout to ensure the appropriate intensity was maintained. 



Participants undertook a second blood test immediately after the 
exercise intervention and a third and final blood test 60 minutes 
after completing the exercise bout (Figure 1). Participants were 
instructed not to consume caffeine, alcohol or nicotine between the 
second and final blood tests. Blood was collected using 10 ml K2E 
EDTA Vacutainer blood collection tubes with BD Vacutainer 
Eclipse Blood Collection Needles (BD Biosciences). Samples were 
kept on ice until WBC isolation and storage. All testing was 
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Figure 2. Exercise-induced changes in T cell populations. T cell populations were measured at each time point and expressed as a relative 
percentage of CD3+ T cells (n = 22) in both CD4+ T cells (A) and CD8+ T cells (B). Relative changes in CD45RA+ and CD45RO+ phenotypes were 
assessed in CD4+ T cells (C) and CD8+ T cells (D) respectively. Error bars indicate SEM. f indicates P = 0.05, indicates P<0.05, ***indicates P<0.001. 
doi:1 0.1 371 /journal.pone.0092088.g002 
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Figure 3. Differential regulation of selected miRNAs in unsorted WBCs. Relative expression of each target miRNA was assessed at pre- 
exercise, post-exercise and 60 minutes post-exercise (n = 18). Whilst only a strong trend was observed for miR-181b (A), significant changes in 
regulation were observed for miR-186 (B), miR-15a (C), and miR-96 (D). All data is expressed relative to an average of RNU44 and RNU48. Error bars 
indicate SEM. "indicates P<0.05 and ""indicates P<0.01, and """indicates P<0.001. 
doi:1 0.1 371 /journal.pone.0092088.g003 



conducted between 7:30 am to 10:30 am to limit circadian 
influence. The time from blood draw to white blood cell isolation 
never exceeded 90 minutes. 

Preparation of Cells 

Whole blood samples were spun at 1000 xg (2250 rpm) (4°C) 
for 15 minutes to separate the plasma and haematocrit sub- 
fractions. The isolated buffy coats were removed and treated with 
red blood cell lysis buffer, spun at 300 xg for 10 minutes and 
washed twice in sterile phosphate buffered saline (PBS). This 
whole blood separation technique isolated WBCs, a heterogeneous 
mix of neutrophils, basophils, eosinophils, lymphocytes and 
monocytes. The WBCs were re-suspended in Iscove's Modified 
Dulbecco's Medium (Life Technologies) containing 10% Fetal 
Bovine Serum (FBS) (Life Technologies) and 200 ul (10% of total 
end volume) of Dimethyl Sulfoxide (DMSO) (Sigma-Aldrich). 
Each sample was stored at — 80°C for 24 hours before being 
transferred to liquid nitrogen storage. 

Frozen cells were thawed and incubated for 30 minutes at 4°C 
with the following monoclonal antibodies (mAbs): anti-CD3- 
FITC, anti-CD4-V450, anti-CD8-APC, anti-CD45RA-PEcy7 and 
anti-CD45RO-PE (BD Biosciences). The cells were then washed 
twice with PBS with 1% fetal bovine serum (PBS/FBS) and re- 
suspended in 200 nl of PBS/FBS and 5 \ll of Propidium iodide (PI) 
prior to analysis. All mAbs were titrated to determine optimal 
concentration. Pre-acquisition compensation using single stains 
and fluorescence minus one stains was conducted to remove 
overlapping fluorescence emission spectra. 



Flow Cytometry 

Flow cytometry was performed using a FACSARIA II Flow 
Cytometer (BD Biosciences), utilizing a red laser emitting at 633- 
nm, a blue laser emitting at 488-nm and a violet laser emitting at 
405-nm. Total lymphocytes were electronically gated based on 
forward scatter (FSC-A)/side scatter (SSC-A) distribution. Live 
lymphocytes were gated according to their expression of PI against 
FSC-A. The CD3+ T cell population was identified by FSC-A 
against FITC expression and then further separated into CD3+ 
CD4+ and CD3+CD8+ subsets based on positive antigen 
expression. The CD45RA+ (naive) and CD45RO+ (memory) 
subsets were identified based on differential expression of 
CD45RA and CD45RO. Cells exclusively gated on positive 
expression of CD45RA will encompass naive and CD45RA+ 
effector memory phenotypes in both CD4+ and CD8+ subsets. 
Similarly, exclusively gating on positive expression of CD45RO 
will encompass central memory and effector memory phenotypes 
[31]. Sorted cells were sorted into PBS and stored at — 80°C. All 
data were processed using Flowjo flow cytometry analysis software 
(Tree Star). 

RNA Extraction 

Total RNA was extracted using TRIzol (Life Technologies) 
according manufacturer's instructions. All RNA samples were 
quantified by spectrophotometry using a Nanodrop (Thermo 
Fisher). 
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Figure 4. Differential regulation of selected miRNAs in sorted T cell subset pools, (n = 22): Each miRNA was assessed in T cell subset pools. 
miR-181b was expressed in CD4+CD45RA+ and CD4+CD45RO+ T cells (A) and CD8+CD45RA+ and CD8+CD45RO+ T cells (B). miR-186 was also 
expressed in CD4+CD45RA+ and CD45RO+ T cells (C) and CD8+CD45RA+ and CD8+CD45RO+ T cells (D). 
doi:1 0.1 371 /journal.pone.0092088.g004 



miRNA Expression Microarrays 

Genome-wide miRNA expression arrays (Agilent Human 
miRNA Microarray, Release 19.0) were performed on pre- 
exercise and immediately post-exercise samples from a subset of 
10 male participants closely matched for age, BMI and V0 2pea j £ . 
The microarrays were performed at the Ramaciotti Centre for 
Gene Function Analysis (University of New South Wales, Sydney, 
Australia), as previously described [32] . The data set obtained has 
been deposited in the NCBI Gene Expression Omnibus database 
according to the Minimum Information About a Microarray Experiment 
(MIAME) guidelines [33], with series accession number 
GSE45041. 

Selection and Validation of Candidate miRNAs 

Differentially expressed miRNAs were further analysed for 
potential binding to telomeric gene transcripts using miRGEN 
Targets (http:/ / www.diana.pcbi.upenn.edu/ cgi-bin/ miRGen/ v3/ 
Targets.cgi) [34], which simultaneously collates and analyses 
intersections and unions of prominent in silico analyses such as 
PicTar, miRagen, microRNA.org, miRanda and TargetScan. 
Four miRNAs: miR-181b, miR-186, miR-15a and miR-96 were 
selected for individual validations via qPCR based on predicted 
interactions between the following miRNAs and telomeric gene 
transcripts: miR-181b and TERT, miR-186 and TERF2IP, RAD50 
and SIRT6, miR-96 and TERF2IP and miR-15a and TATA box 
binding protein (TBP) (Table SI). 

TaqMan assays (Life Technologies) were used to validate the 
target miRNAs in 18 male participants at the pre- and post- 
exercise time-points in accordance with the miRNA expression 



arrays (Table S2). Additionally the target miRNAs were assessed in 
the same 1 8 males at the 60 min post-exercise time-point. Briefly, 
250 ng of total RNA was reverse transcribed for primers using 
TaqMan MicroRNA Reverse Transcription kit (Life Technolo- 
gies) according to manufacturer's instructions. All reactions were 
performed in a BioRad thermocycler (BioRad). 

The qPCR reactions were performed in duplicate in a Viia7 
Real-Time PCR System (Applied Biosystems). All reactions were 
normalized to the average of RNU44 and RNU48; both of which 
have been used extensively as endogenous controls in exercise and 
immunological studies [21,22,35]. Validations were also per- 
formed in pooled CD4+CD45RA+ T cells, CD4+CD45RO+ T 
cells, GD8+CD45RA+ T cells and CD8+CD45RO+ T cells using 
the above method. 

Global Expression of Telomere Extension Genes 

To assess the acute effects of exercise on a wide range of 
telomeric genes, pooled WBC RNA from each time point was 
analysed using a TaqMan Array Human - Telomere Extension by 
Telomerase (Life Technologies). Each plate contained 28 assays 
specific to telomere extension by telomerase associated genes and 
four assays to candidate endogenous control genes; all reactions 
were performed in triplicate. Genes were selected for individual 
sample validation based on fold difference between the three time- 
points. 

Validation of Candidate Telomeric Genes 

Candidate telomeric genes were assessed at three time points via 
qPCR. Total RNA was reverse transcribed using the Applied 
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Figure 5. Differential regulation of 7£7?7"mRNA expression. TERT mRNA expression was assessed in unsorted WBCs (n = 17) (A), CD4+CD45RA+ 
and CD4+CD45RO+ T cells (pool of n = 22) (B), and CD8+CD45RA+ T cells (pool of n = 22) (C). TERT was not detectable in CD8+CD45RO+ T cells. Gene 
expression data is expressed relative to endogenous reference gene (GAPDH). Indicates P<0.05, **indicates P<0.01. 
doi:1 0.1 371 /journal.pone.0092088.g005 



Biosystems High Capacity Reverse Transcription Kit (Life 
Technologies). The qPCR reactions were performed for TERT, 
sirtuin 6 (SIRT6), RAD50 homolog (S. cerevisiae) (RAD50) and 
telomeric repeat binding factor 2, interacting protein (TERF2IP, 
also known as RAP1) in a Viia7 PGR System (Life Technologies). 
Details of qPCR primers are listed in Tables S3 and S4. Target 
genes were normalized to glyceraldehyde-3-phosphate dehydro- 
genase [GAPDH), and analysed using the 2" AA Ct method [36]. 

Pooled T Cell Subsets 

To identify transcriptionally responsive subsets, CD45RA+ and 
CD45RO+ subsets in both CD4+ and CD8+ T cells were assessed 
at each time point. A total of 1 00 ng of total RNA was pooled 
from each cell population at each time point from 22 male 
participants. Sample pooling identifies the transcriptional charac- 
teristics of specific cell populations as opposed to individuals and 
reduces the effects of biological variation. The small RNA yield 
from the sorted cell subpopulations precluded large numbers of 
individual samples. Whilst this technique restricts the scope of 
stringent statistical analysis, it does provide a transcriptional profile 
that can be compared to unsorted WBCs. Reverse transcription 
and qPCR reactions were performed according to previously 
outlined protocols. 

Statistical Analysis 

miRNA microarray samples were between-array normalized 
using the quantile method in Partek Genomics Suite (version 6.6). 
Differentially expressed miRNAs were identified using a paired t- 
test false discovery rate (FDR) <0.05. qPCR data were assessed 



using Friedman's repeated measures for non-parametric data and 
repeated measures ANOVA with post-hoc for parametric data. 
Statistical significance was set at P<0.05. All statistical analysis was 
performed using SPSS Version 17. 

Results 

Exercise Testing and Intervention 

The treadmill test duration was 12. 5 ±1.6 minis. The 0 2 peak 
achieved was 49.3±4.7 mL-kg-min corresponding to the 75 th 
percentile of maximal aerobic power for males aged 20-29 years. 
The maximal respiratory exchange ratio (RER) achieved during 
the treadmill test was 1.2±0.1. The work rate performed during 
the the 30 minutes of treadmill running corresponded to 
80.8±7.1% of O z peak and the average RER was 1.0±0.1. A 
summary of the treadmill ramp test and 30 minute exercise 
intervention data appears in Table 2. 

Lymphocyte Response to Exercise 

Relative frequencies of CD4+ T cells (expressed as a percentage 
of CD3+ T cells) decreased from 43.7% pre-exercise to 36.7% 
post-exercise (P<0.001) and increased to 48.2% 60 min post- 
exercise (P<0.001) (Figure 2). The relative percentage of CD8+ T 
cells (expressed as a percentage of CD 3+ T cells) underwent a non- 
significant increase from 41.9% pre-exercise to 44.8% post- 
exercise before significantly decreasing to 39.9% 60 min post- 
exercise (P<0.01). The relative frequency of CD4+CD45RA+ T 
cells (expressed as a percentage of CD4+ T cells) decreased from 
49.5% post-exercise to 46.8% 60 min post-exercise (P = 0.05). The 



PLOS ONE | www.plosone.org 



6 



April 2014 | Volume 9 | Issue 4 | e92088 



Exercise and Telomere Homeostasis 



□ CD4+CD45RA+ T cells DCD4+CD45RO+ T cells 



3.50 
3.00 
2.50 
2.00 
1.50 
1.00 
0.50 
0.00 



12 
10 
8 

e 

4 

2 
0 



Pre-exercise 



Post-exercise 60min post-exercise 
Time point 



□ CD8+CD45RA+ T cells HCD8+CD45RO+ T cells 



Pre-exercise Post-exercise 60m in post-exercise 

Time point 



£ 



to 



18 
16 
14 
12 
10 
8 
6 
4 
2 
0 



Pre-exercise Post-exercise 60min post-exercise 

Time point 



Figure 6. Differential regulation of 5//?7SmRN A expression. In Unsorted WBCs (n = 17) (A), CD4+CD45RA+ and CD4+CD45RO+ T cells (pool of 
n = 22) (B), and CD8+CD45RA+ and CD8+CD45RO+ T cells (pool of n = 22) (C). Gene expression data is expressed relative to endogenous reference 
gene (GAPDH). "indicates P<0.05. 
doi:1 0.1 371 /journal.pone.0092088.g006 



relative frequency of CD8+CD45RA+ T cells decreased from 
56.3% pre-exercise to 52.1% post-exercise (P<0.05) before 
returning to pre-exercise levels 60 min post-exercise. There were 
no significant changes in either of the CD45RO+ subsets. 

The Acute Effect of Cardiorespiratory Exercise on 
Genome Wide miRNA Expression 

Fifty-six miRNAs were significantly differentially regulated in 
ten healthy males after 30 minutes of intense cardiorespiratory 
exercise (Table S5). 

qPCR Validations of Selected miRNAs 

The qPCR validations (n = 18) identified a non-significant, post- 
exercise upregulation trend in miR-15a (fold change = 1.25), miR- 
181b (fold change = 1.49), miR-186 (fold change = 1.18) and miR- 
96 (fold change = 1.14) (Figure 3). Statistically significant upregu- 
lation between pre- and 60 min post-exercise was observed in 
miR-186 (fold change =1.93, P<0.001), miR-15a (fold 
change = 3.61, P<0.001) and miR-96 (fold change = 2.63, P< 
0.001) (Figure 3). Statistically significant upregulation also 
occurred between post- and 60 min post-exercise in miR-186 
(fold change =1.58, P<0.01), miR-15a (fold change = 3.04, P< 
0.001), and miR-96 (fold change = 2.11, P<0.01). 

The Effect of Exercise on T Cell Subset miRNA Expression 

Only miR-181b and miR-186 were detected in the sorted T cell 
subsets. The expression profile of miR-181b demonstrated a 
biphasic post-exercise downregulation in CD4+CD45RA+ and a 
marginal post-exercise increase in CD8+CD45RA+ T cells 



(Figure 4). There was no appreciable regulation in miR-181b for 
CD4+CD45RO+ T cells whilst CD8+CD45RO+ T cells exhibited 
a stepwise increase. Expression of miR-186 exhibited a 60 min 
post-exercise increase in CD4+CD45RA+ T cells and CD8+ 
CD45RO+ T cells showed a post-exercise stepwise increase. 

Telomere Gene Pathway Expression in Pooled Samples 

Fold changes in relative expression appear in Table S6. The 
TaqMan Human Telomere Extension array was only performed 
once due to limited sample; therefore results are limited to fold 
change without corresponding measures of statistical significance. 
In response to acute exercise, 1 6 of the 28 telomeric genes were 
upregulated from pre- to post-exercise and 15 were upregulated 
from pre- to 60 min post-exercise. A total of 10 genes were 
marginally downregulated between pre-post exercise and 1 3 were 
downregulated from pre- to 60 min post-exercise. 

Validation of Gene Expression in Individual Subjects 

We further validated expression levels of TERT mRNA by 
qPCR based on the fold change and the important role it plays in 
telomere homeostasis. We also investigated the expression of 
SIRT6 mRNA (not featured on the TaqMan Human Telomere 
Extension array) as this gene plays a role in telomeric chromatin 
maintenance [37]. Despite low basal expression levels, WBC 
TEitTmRNA expression exhibited a significant increase from pre- 
to 60 min post-exercise (fold change = 19.4, P= 0.001) (Figure 5). 
To determine if a particular T cell subset was driving the observed 
changes, population specific pools were assessed over the three 
time points. The stepwise upregulation trend was broadly 
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confirmed in CD4+CD45RA+, CD4+CD45RO+ and CD8+ 
CD45RA+ T cell subsets (Figure 5). 

SIRT6 mRNA expression showed significant upregulation in 
WBCs between pre- and 60 min post-exercise (fold change = 1.67, 
P<0.05) and between post- and 60 min post-exercise (fold 
change =1.66, P<0.05). SIRT6 was downregulated immediately 
post-exercise in CD4+CD45RA+, CD8+CD45RA+ and CD4+ 
CD45RO+ T cell subset pools (Figure 6). 

miRNA Target Gene Prediction 

We used miRGen Targets prediction software to determine 
potential target gene transcripts of the selected miRNAs. We 
chose to validate the expression profiles of the potential miRNA 
target gene transcripts for TERF2IP and RAD50 using qPCR. The 
pre- to 60 min post-exercise upregulation of the selected miRNAs 
did not result in statistically significant downregulation of the 
transcript targets over the same time period. The downregulation 
of miR-181b expression between post- and 60 min post-exercise 
showed marginal significance (fold change = 0.84, P = 0.05) and 
paralleled significant upregulation of potential target TERT 
transcript (fold change = 14.75, P<0.05). Expression increases in 
miR-186 from post- to 60 min post-exercise (fold change = 1.58, 
P<0.001) paralleled a significant upregulation of SIRT6 mRNA 
(fold change =1.66, P<0.05), a significant downregulation of 
TERF2IP mRNA (fold change = 0.54, P<0.01) and a strong trend 
towards significance for RAD50 (fold change = 0.76, P=0.05). 
Expression increases in miR-96 from post- to 60 min post-exercise 
(fold change = 2.11, P<0.01) also paralleled simultaneous down- 



regulation in TERF2IP mRNA expression (fold change = 0.54, P< 
0.01) (Table S7). 

White blood cell expression of TERF2IP mRNA was upregu- 
lated immediately post exercise (fold change = 1.46, P<0.01); but 
showed significant downregulation from post- to 60 min post- 
exercise (fold change = 0.54, P<0.01) (Figure 7). TERF2IP mRNA 
expression demonstrated a stepwise downregulation from pre- to 
60 min post-exercise in CD4+CD45RA+ and CD8+CD45RA+ T 
cell subset pools. RAD50 exhibited a non-significant post-exercise 
upregulation of mRNA expression in WBCs followed by a return 
to resting levels (P = 0.28). A stepwise downregulation from pre- to 
60 min post-exercise was observed in CD4+CD45RA+ and CD8+ 
CD45RA+ T cell subset pools (Figure 8). 

Discussion 

Here we report for the first time that acute exercise can lead to 
the transcriptional regulation of several key telomeric genes in 
immune cells. First, we report the upregulation of TERT mRNA 
which plays a critical molecular role in telomere maintenance 
[38]. Second, we report that exercise regulates miRNAs with the 
potential to control the downstream expression of genes involved 
in telomere homeostasis. The exercise-induced regulation of 
telomeric genes and miRNAs may provide an important 
mechanistic link between physical activity, telomeres and im- 
proved health. This is an important finding given the extensive 
evidence linking accelerated telomere shortening to several 
chronic diseases [39], many of which can be ameliorated by 
aerobic exercise [40]. Telomere homeostasis underpins the 
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function of several immune cell subsets [41], which in turn play 
critical roles in chronic pathology such as age-related diseases [42], 
atherosclerosis [43] and metabolic diseases [44]. 

Exercise-induced miRNA regulation exhibits exquisite specific- 
ity for both tissue type and exercise modality. Differential 
expression profiles have been identified between muscle contrac- 
tion types [45], aerobic fitness levels [23] and resistance training 
adaptability [46] . Recent studies have assessed miRNA expression 
in leukocyte subsets, identifying differential regulation of 38 
miRNAs in neutrophils [22], 23 miRNAs in natural killer cells 
[35] and 34 miRNAs in peripheral blood mononucleated cells 
(PBMCs) [21,22,35]. Thirty eight miRNAs were differentially 
expressed in neutrophils after ten 2 minute bouts of intense cycle 
ergometry [22]. Pathway analysis of miRNA targeted genes 
identified pathways responsible for neutrophil immune function 
and apoptosis [47,48], chronic inflammation [49] and immune 
function [50]. The analysis of PBMC miRNA revealed 34 
significantly regulated miRNAs that influenced genes associated 
with 1 2 signalling pathways including pro- and anti-inflammatory 
cytokine regulation [51-53], lymphocyte activation and differen- 
tiation [54], cell communication [55] and cancer [56]. 

The current study provides a unique and unprecedented snap 
shot of over 1 300 miRNAs, revealing 56 significantly regulated in 
response to exercise. In silico analyses predicted miRNA/mRNA 
interactions between the following miRNAs and telomeric gene 
transcripts: miR-181b and TERT, miR186 and TERF2IP, RAD50 
and SIRT6, miR-96 and TERF2 IP and miR-15a and TATA box 
binding protein (TBP). Validations conducted for TERT and 
SIRT6 showed significant 60 min post-exercise increases that 



paralleled concomitant increases in the potential binding miRNA, 
effectively dismissing any significant miRNA/ transcript interac- 
tion. Validations for TERF2IP and RAD50 revealed significantly 
decreased transcript abundance at 60 min post-exercise for 
TERF2IP only, paralleled by concomitant increases in potential 
binding miRNAs (miR-186 and miR-96). 

Two miRNAs investigated in the present study; miR-181b and 
miR-96, have previously been identified in neutrophils following 
30 minutes of exercise [22]. Expression profiles of miR-181b were 
broadly confirmed at the same time point between the two studies; 
however, miR-96 was downregulated in the previous study and 
upregulated in the present study. This may be due to analysis of 
whole WBCs in the present study (including neutrophils) as 
opposed to isolated neutrophils in the previous study. Pathway 
analysis identified 578 genes targeted by miR-96 and 652 genes 
targeted by miR-181b [22]. A more recent study identified 
upregulation of miR-181b and miR-15a in PBMCs immediately 
after 30 minutes of intermittent aerobic exercise [21]. Both miR- 
181b and miR-15a were identified in the present study; however, 
only miR-15a was significantly regulated. Utilizing a similar 
exercise protocol, a recent study identified four miRNA-mRNA 
networks dynamically regulated by 30 minutes of exercise [57]. 
The target mRNAs are involved in apoptosis, immune function, 
transcription regulation and membrane traffic of proteins. 

RAD50 associates with meiotic recombination 1 1 (MRE1 1) and 
nibrin (NBS1) to form the MRN complex which repairs DNA 
damage and assists telomere structure [58,59]. The MRN complex 
also positively regulates telomerase-dependent telomere elongation 
via an interaction with telomeric repeat binding factor 1 (TERF1) 
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and ataxia telangiectasia mutated (ATM) [60]. The complex 
formed by RAD50 and MRE1 1 is thought to help stabilize t-loop 
formation [61]. There is no immediately intuitive reason for the 
post-exercise downward trend in RAD50 mRNA expression. 
Possible hypotheses include transiently compromised stability of 
the shelterin complex and suppression of DNA damage machinery 
at the telomere by the heavily fortified shelterin complex [62]. 
Telomeric DNA is preferentially damaged by oxidative stress [63]; 
however, the conformation of the shelterin complex and T-loop 
may preclude access to sites of DNA damage. Mass spectrometry 
and affinity purification have also identified an association 
between RAD50 and the TERF2 IP /telomere repeat binding 
factor 2 (TRF2) protein complex [64]. This association suggests 
that RAD50 mRNA expression may parallel that of TERF2IP 
mRNA. 

The hypothesis of exercise-induced shelterin instability is 
reinforced by the simultaneous decrease in TERF2IP mRNA, 
the protein of which forms part of the shelterin complex and is 
recruited to telomeres via interaction with TRF2 (also known as 
TERF2) [65]. TERF2IP deficiency reduces telomere stability and 
increases recombination [66]. Equivocal results have identified 
TERF2IP as both a negative regulator [64] and a positive 
regulator of telomere length [65] . Additional roles for TERF2IP 
include prevention of non-homologous end joining [67] and 
homology-directed repair [66], protection from obesity via 
regulation of metabolic genes [68] and regulation of senescence 
[69]. Additionally, TERF2IP associates with iKB kinases in the 
cytoplasm and regulates NF-kB modulated gene expression [70]. 

To our knowledge, this study is the first to report the 
upregulation of TERT mRNA after acute aerobic exercise in 
WBCs and T cell subsets. The current findings may provide a 
potential mechanistic link between physical activity and telomere 
length. Transcriptional plasticity of TERT is important for the 
adaptive immune system; the replicative capacity of which 
depends upon the telomere length and TERT expression of the 
constituent T cells [71]. Ectopic expression of TERT has been 
shown to increase CD4+ and CD8+ proliferative capacity and 
enhance resistance to oxidative stress and apoptosis [72-74]. 
Unlike T cells, mature granulocytes of the innate immune system 
do not undergo successive rounds of division therefore granulocyte 
telomere length is a function of myeloid progenitor cells [75]. 
Accordingly, telomerase activity within the granulocyte population 
is low [76] to undetectable [77,78]. Whilst the WBC fraction 
analysed in the present study contained elements of both the 
adaptive and innate immune system, all findings were subsequent- 
ly validated in sorted populations of T cells subsets. 

Epel et al., (2010) identified an 18% increase in leukocyte 
telomerase within one hour of an acute psychological stressor [79]. 
The increase was independent of immune cell redistribution and 
was positively associated with concomitant increases in Cortisol. 
The authors postulated a possible protective function for the acute 
increase in telomerase, such as preparation for immune cell 
proliferation [79]. The authors also indicated that phosphorylation 
of telomerase and/or changes in subcellular localization likely 
account for the acute increases in telomerase, given the extended 
timeline likely required for transcriptional regulation and/or 
alternate splicing of TERT [79] . The findings of the present study 
are in opposition to this hypothesis, demonstrating transcriptional 
regulation of TERT within 60 minutes of exercise. 

There is significant evidence of a linear relationship between 
TERT mRNA expression and telomerase activity in a range of 
tissues [80-84]. Despite this, numerous human tissues exhibit 
significant heterogeneity between the two, indicating potential 
post-transcriptional modification of TERT. Gizard et al, (2011) 



found that macrophages stimulated with lipopolysaccharide, 
oxidized low density lipoprotein and tumour necrosis factor-a 
exhibited induction of TERT mRNA expression that paralleled 
inducible telomerase activity [80]. However, macrophages stim- 
ulated with interleukin- 1 (3 exhibited only modest increases in 
TERT mRNA despite high induction of telomerase activity [80]. 
Although demonstrating an overall correlation between TERT 
mRNA and telomerase activity in lung cancer tissues, Hara et al., 
(2001) reported that 12.9% of samples exhibited TERT mRNA 
without telomerase activity and an additional 12.9% demonstrated 
telomerase activity without detectable TERT mRNA expression 
[81]. 

Human thymocytes, tonsil and peripheral blood T and B cells 
have all been shown to exhibit TERT mRNA expression 
independent of telomerase activity [85]. All lymphocyte subsets 
harvested from peripheral blood and thymus contain TERT 
mRNA and protein irrespective of telomerase activity [86]. 
Significant discrepancies between TERT mRNA expression and 
telomerase activity have also been identified in soft tissue sarcomas 
[87] and non-small cell lung cancer [88]. Additionally, expression 
levels of TERT mRNA do not regulate telomerase activity in 
telomerase-reconstituted primary human fibroblast clones [89] . 

The numerous associations between exercise, telomere length 
and telomerase are observational in design with either confound- 
ing lifestyle factors or discrepant blood collection timelines. A 
recent study identified no change in PBMC TERT mRNA 
expression or telomerase levels after 7 days of ultra-marathon 
running in trained athletes [90]. Given that measurement 
occurred the day after marathon completion, acute changes in 
gene expression or T cell frequencies may have regressed to basal 
levels within this time frame. The current study is unique in that it 
used an acute measurement timeline; providing a snapshot of 
transcriptional changes immediately after and 60 min after 
exercise. Parallel changes in cell population frequencies were also 
assessed to account for subset redistribution. 

The acute inflammatory milieu induced by exercise may in part 
be responsible for signalling increases in WBC TERT mRNA. 
Previous research has shown pro-inflammatory signalling induces 
TERT mRNA expression and telomerase activity in macrophages 
[80]. Intense exercise transiendy and preferentially redistributes T 
cell subsets [91,92]. Studies have equivocally identified preferential 
mobilization of T cells with shortened telomeres [93], senescent 
phenotypes with longer telomeres [94] and CD45RA+ effector- 
memory phenotypes [95]. We demonstrated a non-significant 
60 min post-exercise decline in CD4+CD45RA+ T cells and a 
non-significant increase in CD4+CD45RO+ T cell frequency. 
Telomere length was not measured in the present study because 
the focus was acute regulation of telomeric genes. Given the 
average trajectory of leukocyte telomere shortening [96-98], 
genuine changes in telomere length would not be detectable after 
60 minutes. Additionally, telomere length appears sensitive to the 
stage of T cell differentiation [94]; a factor our flow cytometry 
staining panel could not accurately discriminate. 

Exhaustive exercise has been shown to increase human PBMC 
expression of SIRT1 and decrease SIRT3 and SIRT4 mRNA [99]. 
To our knowledge, this is the first study to characterize exercise- 
induced SIRT6 mRNA expression in humans. SIRT6 is a 
chromatin binding protein existing predominantly in the nucleus 
[100]. SIRT6 influences telomeric chromatin by deacetylating 
H3K9 and H3K56 [37], thereby reducing chromatin accessibility. 
SIRT6 also plays an important role in DNA repair mechanisms by 
modulating base excision repair (BER) [101] and double strand 
break (DSB) repair [102]. SIRT6 is recruited to DSB sites under 
conditions of oxidative stress [103]. Other extra-telomeric roles for 
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SIRT6 include enhanced mitochondrial respiration [104], tran- 
scriptional regulation of gene expression [105], positive regulation 
of pro-inflammatory cytokines [106] and systemic glucose 
metabolism [107]. The post-exercise upregulation of SIRT6 
mRNA exhibited in unsorted WBCs was not replicated in sorted 
T cell populations. A likely explanation for the disparate results is 
that the predominant cell type(s) driving the expression changes 
were not T cells and were therefore excluded when T cells were 
positively selected. 

There are some inherent limitations to this study. Without a 
priori knowledge of the precise time course of miRNA expression, it 
is possible that the time course used in this study may have missed 
the greatest magnitude of change. Additionally, the magnitude of 
miRNA expression needed to elicit detectable decreases in target 
mRNA is unknown. Determining the optimal timeframe in which 
to measure telomeric gene expression also presents a challenge as 
littie is known about the exercise-induced transcriptional timeline 
or half-life of telomeric gene transcripts. This study assessed 
changes in mRNA expression and did not quantify protein or 
functional enzyme status. Additionally, we were not able to 
differentiate between the telomeric and extra-telomeric roles of the 
genes and miRNAs in this study. Observed increases in SIRT6 and 
TERT mRNAs may have resulted from the upregulation of extra- 
telomeric pathways. Additional reporter assay and gain/loss of 
function experiments are needed to confirm the legitimacy of the 
miRNA/ mRNA interactions. 

The low sample yield from sorted T cell populations necessi- 
tated pooling into cell/ time point specific pools for analysis. Whilst 
this provided interesting, subpopulation-wide overviews, it pre- 
cluded additional individual validations and robust statistical 
analysis. Whilst the analysis was predominantly conducted in 
WBCs and T cell subsets, analysis of other leukocyte subsets such 
as B cells and neutrophils may help form a more complete picture. 

Cytomegalovirus (CMV) is a persistent and ubiquitous herpes 
virus that can affect T cell telomere/telomerase homeostasis [108- 
110]. The CMV seroprevalance rate in Australian males aged 20- 
24 years is 50% [1 1 1], with global seroprevalence rates exceeding 
70% above 60 years of age [112]. Financial and logistical 
constraints precluded the individual testing for CMV seropositivity 
and we felt that testing and excluding positive individuals would 
place further constraints on an already small sample size. Potential 
participants were excluded if they had been diagnosed with 
glandular fever, chronic fatigue syndrome or other viral infections 
lasting longer than 3 months. Given the asymptomatic nature of 
CMV infection, we are unable to exclude the possibility that some 
participants were unknowingly seropositive. 

In conclusion we have shown that 30 minutes of cardiorespi- 
ratory exercise is sufficient to elicit an upregulation of key 
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